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Abstract: The main aim of this paper is twofold: (1) reveahng a relation between the 
counting function N (A) (the number of the eigenstates with eigenvalue smaller than a given 
number) and the heat kernel K (t), which is still an open problem in mathematics, and 
(2) introducing an approach for the calculation of N (X), for there is no effective method 
for calculating N (A) beyond leading order. We suggest a new expression of N (A) which 
is more suitable for practical calculations. A renormalization procedure is constructed for 
removing the divergences which appear when obtaining (A) from a nonuniformly conver- 
gent expansion of K (t). We calculate A (A) for Z?-dimensional boxes, three-dimensional 
balls, and two-dimensional multiply-connected irregular regions. By the Gauss-Bonnet 
theorem, we generalize the simply-connected heat kernel to the multiply-connected case; 
this result proves Kac's conjecture on the two-dimensional multiply-connected heat ker- 
nel. The approaches for calculating eigenvalue spectra and state densities from A (A) are 
introduced. 
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1. Introduction 

A problem stemming from physics and soon becoming an important mathematical problem 
is to recover geometry of a manifold from the knowledge of the eigenvalues of a natural 
differential operator. This problem originates in the theory of radiation (how to determine 
the state density of the electromagnetic wave in a given cavity) and is formulated by Kac 
as "Can one hear the shape of a drum?" Q. 

What the original problem asked is that how many eigenvalues are smaller than a given 
number A. This problem is formulated as to seek the so-called counting function N (X). 
For a given spectrum {A„}, the counting function is directly defined to be 

(A) = the number of eigenvalues (with multiplicity) smaller than A, (1-1) 
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i.e., 

NiX)=^l. (1.2) 

A„<A 

In other words, the counting function represents the number of the eigenstates whose 
eigenvalues are smaher than A. Nevertheless, the counting function N (A) is very difficult 
to calculate and there is no general method for calculating N (X) in mathematics Q. 
This is because when calculating A^(A), one often encounters some unsolved problems in 
number theory. For example, when calculating the counting function for the spectrum of 
the Laplace operator on a tori, one encounters the Gauss circle problem in number theory. 
The Hardy-Littlewood-Karamata Tauberian theorem gives the first term of the asymptotic 
expansion of (A), but does not provide any information beyond the first-order term H]. 
For a given spectrum {An}, one also introduces the heat kernel, 

Kit) = Y,e-^-\ (1.3) 

n 

which is an another important function describing the relation between the eigenvalues of 
the operator and the geometrical property of the manifold. The heat kernel is relatively 
easy to calculate and some methods for calculating heat kernels have been developed . In 
special, Elizalde et al. developed a very effective approach for calculating K (t) Q |5|, |6|, |^, |^ . 

There is no doubt that there must exist a relation between the counting function N (A) 
and the heat kernel K{t), since they are both determined by the same spectrum {A„}. 
Nevertheless, the relation between (A) and K (t) is still an open problem in mathematics 

Though the definition of counting function (A), eq. (|1.1| ) or ( |1.2[) , is comprehensive 
and intuitive, other than the definition of heat kernel K (t) , eq. ( |1.3| ) , the definition of A'^ ( A) 
is not suitable for practical calculations. In fact, if starting from this definition to calculate 
the counting function, one may always encounter the mathematical difficulty mentioned 
above. In this paper, we will propose an expression for A^(A), which has a similar form 
with the definition formula of K (t), eq. ( |1.3D , and is more operable than the definition 
( |1.2| ). This expression can be regarded as an alternative definition formula for A^ (A). 

As a main result of this paper, we provide a relation between the counting function 
A^ (A) and the heat kernel K (t). This problem interests us from both purely mathematical 
viewpoint and practical viewpoint. This relation allows us to calculate A^ (A) from K (t) 
which, as mentioned above, is often relatively easier to calculate. 

The relation between the counting function A^ (A) and the heat kernel K (t) presented 
in this paper is an integral transformation. However, in practice the heat kernel is often 
given in the form of a series expansion and in most cases only the first several expansion 
coefficients can be obtained This requires us to integrate term by term. Nevertheless, 
the series expansion of K (t) is not convergent uniformly. As a result, the integral of some 
terms will be divergent. The counting function A^ (A) is of course finite; the divergences are 
caused by illegally integrating a nonuniformly convergent series term by term. For dealing 
with this problem, we provide a renormalization procedure to remove the divergences. 



- 2 - 



In fact, the results of this paper provide two approaches for calculating the counting 
function N (X): (1) by the expression for N (X) presented in this paper, one can calculate 
N (X) directly, and (2) by the relation between N (X) and K (t), one can first calculate 
the heat kernel K (t) and then calculate (A) from K (t). As a comparison, we calculate 
N (A) for a D-dimensional box by both these two approaches, respectively. 

As an example of the calculation of counting functions by the relation between N (A) 
and K (t), we calculate N (A) from K (t) for three-dimensional balls. 

As another example, we calculate the counting function from the heat kernel for the 
minus Laplace operator in a two-dimensional region with irregular shape and nontrivial 
topology. For two-dimensional heat kernels, there are two known results: the heat kernel 
for a simply-connected region and a hypothesis made by Kac on the heat kernel for the 
multiply-connected case Q. In the present paper, we first generalize the heat kernel for 
a simply-connected region given in to the case of a multiply-connected region which is 
bounded by a smooth but irregular curve and with some holes. This result proves Kac's 
conjecture. Then, we calculate the counting function for the two-dimensional multiply- 
connected region from the heat kernel by the relation between N (A) and K (t). 

As applications, we discuss the problem that how to calculate the asymptotic expres- 
sions for eigenvalue spectra and state densities from the counting function N{X). By 
(A), we construct an equation for eigenvalues and then provide an approximate solution 
in which the eigenvalue is expressed as a function of heat kernel coefficients. Moreover, we 
also provide two expressions for state densities by the relation between A^ (A) and K (t) 
and by the expression of A^ (A) provided in the present paper. 

The counting function and the heat kernel are interesting in both mathematics and 
physics: 

In mathematics, the relation between the spectrum of the Laplace operator on a Rie- 
mannian manifold and the geometry of this Riemannian manifold is an important subject 
||2|, |lO|, and the problem of spectral asymptotics is one of the central problems 

in the theory of partial diff'erential operators ||l3[. One of the main problems is to seek 
the asymptotic expansions of the counting function A^ (A) and the heat kernel K (t). The 
general relation between A^ (A) and K (t) is still unknown. Especially, there is no general 
method for calculating the asymptotic expansion of the counting function A^ (A). 

In physics, the spectrum of the Laplace operator on a Riemannian manifold can be 
directly applied to boundary quantum fields. Moreover, reconstructing the geometrical 
property of a system from an eigenproblem is an interesting and important problem. For 
example, Aurich et al. reconstruct the shape of the universe from the result of the eigen- 
problem [0]. The counting function is directly related to the Casimir effect, the spectrum 
problem, and the state density, etc. Moreover, there are many studies on heat kernels. 
In quantum field theory, it is of crucial importance to evaluate the one-loop divergences, 
and the control of the ultra-violet divergences can be achieved by using heat kernel reg- 
ularization methods HI, [l^, 16]. The heat kernel expansion becomes a standard tool in 
the calculations of vacuum energies |17|, the Casimir effect [18, and quantum anoma- 
lies |2C]. Now a large amount of research has been devoted to quantum gravity based on 
the heat kernel expansion, including semiclassical approaches |21], black hole thermody- 
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namics [22|, wormhole physics [23|, heat kernels in curved space |24|, and supergravity 
theories [25, 26]. In addition, the heat kernel expansion is also used to study string theory 
1 27, 28, ^ and noncommutative field theory pO| , 32 1. The relation between the counting 
function and the heat kernel allows one to introduce the concept of the counting function 
into these fields. Moreover, there exist some indirect ways to construct the heat kernel 
without solving the heat equation, and the heat kernel coefficients for various manifolds 
have been discussed |^, ^, 33, 34 1, which makes the heat kernel expansion a very convenient 
and effective tool. By the relation provided in the present paper, such results of the heat 
kernel can be directly applied to the problem of the counting function. 

In section we provide an expression for the counting function (A), which is more 
operable compared to the definition ( |1.2| ) and can be regarded as an alternative defini- 
tion formula for the counting function. In section ^, we point out a relation between the 
counting function N (X) and the heat kernel K(t). In section we provide a renormal- 
ization procedure for removing the divergences appearing in the calculation of N (A) from 
a nonuniformly convergent series expansion of K{t). In section we compare the two 
approaches for calculating N (A) using a Z)-dimensional box as an example. In section 
|6|, we calculate the counting function of a three-dimensional ball from the heat kernel. 
In section ^, for two-dimensional cases, we first generalize the heat kernel in a simply- 
connected region to a multiply-connected region, and then calculate the counting function 
in a multiply-connected region from the heat kernel. The result of the heat kernel in the 
multiply-connected region obtained in the present paper proves Kac's hypothesis. In sec- 
tions ^ and ^ as applications of the counting function, we discuss how to calculate the 
asymptotic expressions for spectra and state densities from A^(A). The conclusions are 
summarized in section ITO . 



2. An expression for the counting function: an alternative definition for- 
mula 

The definition formula for the counting function N (A), eq. ( |l.2D , is not suitable for practi- 
cal calculations. In fact, starting from this definition to calculate (A), one may encounter 
many difficulties, e.g., the unsolved problems in number theory 0. Contrarily, the defini- 
tion formula for the heat kernel, eq. (|1.3| ), is more operable. In this section, we present an 
expression for N (X), which has a similar form with that of the heat kernel K (t), eq. (|1.3|). 
Such an expression can be regarded as an operable definition formula for N (A). 



Theorem 1 

A^(A) = lim V— — i- . (2.1) 

n 

Proof. Observing that 



lim 



1, when A„ < A, 



/3^oo e/^(^"-^) + 1 10, when A^ > A, 
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we have 

lim y -T7r-K\ — = y 1 = iv (A) . 

n A„<A 

Comparing with eq. (|L^), eq. (l2.l|) is obviously more operable, since eq. ( |2.lD converts 
the partial sum in eq. (|1.2|) , X]a„<A' i'^^o ^ ^^^^ possible values, Yl\„<<x- This 
will, of course, make the calculation easier. 



3. The relation between counting function and heat kernel 

In this section, we provide a relation between the counting function N (A) and the heat 
kernel K (t). This relation is an interesting mathematical result and is useful for practical 
calculations. As mentioned above, the calculation of the counting function N (A) is more 
difficult than that of the heat kernel K (t). With the relation between (A) and K (t) given 
in the following, one can first calculate K (t) and then calculate N (A) from the result of 
Kit). 

The relation between (A) and K (t) is as follows. 
Theorem 2 

poo 

K{t)=t / A^(A)e-^*dA. (3.1) 

JO 

Proof. The generalized Abel partial summation formula reads 

B (n) /' (n) dn, (3.2) 

where Aj G R, Ai < A2 < • • • < An < • • • , and lim^^oo An = 00. / (n) is a continuously dif- 
ferentiable function on [tti, 112] (0 < ui < n2, Ai < U2), h (n) (n = 1, 2, 3, • • • ) are arbitrary 
complex numbers, and B [u) = Y2\„<u^ i''^)- ^^Pply the generalized Abel partial sum- 
mation formula, eq. (|3.2|) , with / (u) = e~"(*~'^o) and b{n) = ane~^"'^'^, where s, sq € C. 
Then 

A {U2, s)-A {ui,s) = A {U2, so) e-"2(^-^o)-^ {ui,so) e""i(^-^o)+(s - sq) A (u, so) e-^'^'-^^Uu, 

J ui 

(3.3) 

where 

A{u,s) = a„e-^"^ (3.4) 

\„<u 

Setting On = 1 in eq. (|3.4|), we find 

A(A,0) = ^ 1 = A^(A), 

A„<A 

the counting function, and 

^(oo,t)=^e-^"* = K(t), 

n 



- 5 - 



the heat kernel. By eq. (p^), we also have A{0,t) = 0. Then, by eq. (|3.3| ), we have 

/•oo 

K{t) = A{oo,t)-A{0,t)=t N{X)e^^^dX. (3.5) 

Jo 

This is just eq. ( |3.1| ). 
Theorem 3 

1 rc+icx> Xt 1 

N{X) = — K{t)—dt, c> Urn 4^. (3.6) 

27rz 7c-ioo i '^n 

Proof. By the Perron formula, we have 

1 pc+ioo t 



where 



oo 



n=l ^" 

and c is a constant which is greater than the abscissa of absolute convergence of the Dirichlet 
series / (s). Setting 

a„ = 1 and /x„ = C^"- 



oo 



in eq. ( |3.8D , we obtain the heat kernel, 

f{t) = Y. e-^"* = K{t). 

n=l 

The abscissa of absolute convergence of / (t) equals its abscissa of convergence, equaling 



lim„^oo lnn/A„ = lim„_^oo lnn/A„. Thus, by eq. (3/7), we have 

fc+ioo 

1 



1 pc+too , 



and c > lim^^oo lnn/A„. This proves the theorem. 

Clearly, the relation between the counting function (A) and the heat kernel K (t) 
obtained above is an integral transformation and the corresponding inverse transformation. 

4. Calculating (A) from an asymptotics of K (t) 

In principle, we can calculate the counting function N (A) from the heat kernel K (t) or 
calculate K (t) from (A) by the relation ( ^.6| ) or ( |3.1[ ) directly. However, the relation 
given above is an integral transformation. Therefore, when the heat kernel K (t) is in 
the form of a series expansion, practically, the above relation is available only when this 
integral transformation can be applied to each term of the series, or, in other words, this 
series must can be integrated term by term. If the series is not uniformly convergent, 
even though the integral of the sum function is convergent, the integrals of some of the 
terms may be divergent. In this case, we need a renormalization procedure to remove the 
divergences. 
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4.1 An asymptotics for N (A) 

Concretely, the expansion of heat kernel K (t) can be expressed as 



K(t) = (47rt)-^/2 V Bkt", (4.1) 



fe=o,§,i,- 



where -B^ is the heat kernel coefficient, D is the dimension of space. In fact, in practice 
only the first several heat kernel coefficients can be obtained, so we have to integrate term 
by term. Substituting eq. ( [4.1[ ) into eq. (|3.6| ) and exchanging the order of integration and 
summation gives 

oo 

iV(A) = (47r)-^/2 ^ BklDik), (4.2) 
fc=o,i,i,- 

where 

1 pc+ioo At 

lD{k) = ^ t-^l^+^^—dt. (4.3) 

Nevertheless, for a power series, unless the series is uniformly convergent, it is not per- 
missible to integrate term by term, i.e., the order of integration and summation cannot 
be exchanged. When calculating an integral of a series that is not uniformly convergent 
term by term, one may find that the integrals of some terms diverge. In our case, it can be 
directly seen that only when k < D/2 + 1, (k) is convergent. That is to say, the term 
with k > D/2 + 1 needs to be renormalized for obtaining a finite result. 

4.2 Renormalization 

The integral (|4.3|) is divergent when k > D/2 + 1. For obtaining the counting function 



(A) from the series (4.2), we need to make sense of these divergent integrals. In other 
words, we need a renormalization procedure to remove the divergences. 

It can be directly seen that eq. (|4.3|) is the inverse Laplace transformation of the 



function ^-^/s+fe-i^ ^hus for k < D/2 + 1, we have ll 

The definition domain for the gamma function F (1 + D/2 — k) is k < D/2 + 1. Now, 
we deal with the case of k > D/2 + 1. When —D/2 + /c — 1 is a non-negative integer, 
the integral (U) equals (J^^-^/^+fc-i) ^)^^ where 5^") (A) is the nth derivative of 6 (A). 
When —D/2 + A; — 1 is a positive half integer, we can redefine the gamma function by 
analytically continuing the gamma function via the recurrence relation F (x) = ^F (x + 1) 
j36U . Concretely, when defining the gamma function as F (x) = ^F {x + 1), the definition 
domain of the gamma function is changed to k < D/2 + 2. Repeating this procedure will 
eventually analytically continue the gamma function to the whole real axis. 

From the above result we can see that, like that in the perturbative expansion of the 
scattering amplitude in quantum field theory, the first several terms are convergent (like 
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tree diagrams) and the other terms are divergent and needed to be renormahzed (hke loop 
diagrams). For clarity, write the counting function N (X) in two parts: 



N (A) = Nt (A) + Nl (A) , (4.5) 

where 



-+- 

2^2 



7VT(A) = (47r)-^/2 Bklnik) (4.6) 

is the convergent part and 



iV^(A) = (4^)-^/2 ^^^D{k) (4.7) 



is the renormahzed divergent part, where {k) denotes the renormahzed result. (A) 
is only a higher-order contribution to the counting function and is often negligible. 
After renormalization, eqs. ( |4.6| ) and ( [4.7D can be written as 

oo oo 

N{X)= C,X^^'-'+ E (4vr)-^/^i?i+^/2+,5(')(A) 

fe=0,i,l,- 1=0,1,2,- 
D 

2 OO CO 

= E CkX''/^-''+ Yl CkX^^^-'+ Y (4vr)-^/2^fc5('=-(T+i))(A) 

fc=Oi,l,- '==f + if + l'- fc=f+l,f+2,- 

(4.8) 

with the counting function coefficient 

C, = (47r)-^/2 Bk_ ^ 

" ^ ^ r(l + Z)/2-A:) ^ ^ 



It should be emphasized that the delta-function terms in eq. ( [4.q ) will contribute the 
zero-point energy; this is, such contributions will play an essentially important role in the 
calculation of the Casimir effect. 



5. Comparison of the two approaches for calculating A^(A): D-dimensional 
boxes 

The results presented in sections ^ and |^ show that there are two approaches for calculating 



the counting function A^(A). The first is to directly calculate N {X) by eq. (2A_), the 
operable expression of N (A) presented in this paper; the second is to first calculate the 
heat kernel K (t) and then to calculate N (A) from K (t) by the relation between N (A) 
and K (t). In this section, as an example, we calculate N (A) by the above two approaches 
respectively for a rectangular L'-dimensional box. 

The spectrum of the minus Laplace operator of a D-dimensional rectangle box reads 
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(^1 ^2 A 
Y2 ~'~ T2 + ■ ■ ■ + T2~ I ' ''^i = li 2 • • • , (5-1) 

where Lj is the side length. 

First, we directly calculate the counting function N {X) from eq. ( |2.1| ): 

iV(A)=lim i . (5.2) 

ni,n2,--- ,™d ^ ^ ' -r ± 

oo 

By the Euler-MacLaurin formula = F{n)dn+\F{{))-^B2F' {{))-^B4^F"' {{))+ 

n=0 

• • • , where B2, B4, ■ ■ ■ are Bernoulli numbers, we achieve 
AT (A) = (47r)-^/2 



r(Z)/2 + l) ^ Li r (Z?/2 + 1/2) 



1/ vr 



(5.3) 



where V = L1L2 is the volume. 

Alternatively, we can first calculate the heat kernel from the definition formula of K (t), 

K{t)= Yl e-^("i'"2'-'"o)*. (5.4) 

ni,n2,--- ,nr) 

A direct calculation by the Euler-MacLaurin formula gives the heat kernel coefficients: 

^ V " V 

Bo = V, B„2 = -V^Y.T.^ ^i=^E 



_^ Li . LiL-i 

D 



1=1 i<j ■' 



B,l2 = {-\r^^l^ J J ■ (5.5) 

By the relation between the counting function coefficient and the heat kernel coefficient, 
eq. ([4.9|) , we also achieve eq. (|5.3|) . 

From the above result, we can see the equivalence between these two approaches. 

6. Calculating A^(A) from K (t): three-dimensional balls 

In this section, as an example, we calculate the counting function N (A) from the series 
expansion of the heat kernel K (t) for a three-dimensional ball. It should be emphasized 
that for calculating the heat kernel K (t), it is not needed to know the eigenvalue spectrum 
in advance p. 

For a three-dimensional ball, from eq. (|4.9D , the expansion coefficients of the counting 
function (A) can be obtained as 
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The heat kernel coefficients for Dirichlet, Neumann, and Robin boundary conditions for 
three-dimensional balls are calculated explicitly in [^j. 

Take the case of Dirichlet boundary condition as an example. Ref. Q gives the first 21 
heat kernel coefficients. From this, we can obtain the first 21 counting function coefficients 
Cjt by eq. (6.1). Only taking the tree-diagram- like part into account, from eq. ( [4. 61 ), we 



have 

N (A) = — i?3A3/2 _ 1 ^RX^/^ - 1 - -1 L-. (6.2) 

^ ^ 97r 4 Svr 48 SlSvr i?AV2 ^ ^ 

It can be directly checked that the tree-diagram-like part provides the main contribution 

to the A^(A). 

7. Multiply-connected cases: two-dimensional K (t) and (A) 

An interesting special case is about two-dimensional counting function and two-dimensional 
heat kernel, which is just the original problem formulated by Kac as "Can one hear the 



shape of a drum?". The first result of this problem is proved by Weyl |37| and then im 



proved by Pleijel |38| and Kac ||l|. In []|], Kac pointed out a special relation between N (A) 
and K (t): in two dimensions, the leading-order term of N (A) is accidentally equal to the 
leading-order term of K (t). However, beyond the leading order, there is no further result. 
Moreover, for the multiply-connected case, Kac made a hypothesis about the topologi- 
cal contribution on the heat kernel K {t): the contribution from the nontrivial topology 
(connectivity) is in proportion to the Euler-Poincare characteristic number. 

In this section, (1) with the help of the relation between N (A) and K (t) given above, we 
give a proof for Kac's result: the leading-order term of N (A) equals the leading-order term 
of K (t); (2) based on the Gauss-Bonnet theorem, we first generalize the simply-connected 
heat kernel K (t) given in to the multiply-connected case (this result proves Kac's conjec- 
ture: in two dimensions, the topology contribution is proportion to the Euler-Poincare char- 
acteristic number), and, then, calculate the two-dimensional N (A) in a multiply-connected 
region with the help of the relation between N (A) and K (t). 

7.1 The leading contribution of (A) and K (t) in two dimensions 

Using the relation between (A) and K (t), we first prove that when the number of eigen- 
states per unit interval (the density of eigenstates) p (A) is a constant. A'' (A) equals K (t) 
in the limit A ^ oo or i — > 0, i.e., in the limit A — > oo or t ^ 0, 

N{X) = k(\) or N(l)=K{t). (7.1) 



The proof is straightforward. In the limit A — > oo or 4 — > 0, the summations can be 
converted into integrals: 

NiX)=^l = (\{\')dX\ (7.2) 

/>oo 

(i) = V e-^"* = / p (A') e-^'*(iA'. (7.3) 

-'0 
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Figure 1: Converting a multiply-connected region to a simply-connected region. 

If p (A) = C, where C is a constant, then 

N{X) = C\, (7.4) 

Kit) = j. (7.5) 

This proves eq. ([7?l|). 

In two dimensions, the state density for the eigenstate of the minus Laplace opera- 
tor is a constant, so the leading contribution of the counting function equals the leading 
contribution of the heat kernel. This is just the case appeared in Kac's work [j^. 

7.2 K (t) and N (A) in Multiply-connected regions 

For achieving the counting function in a multiply-connected region by the relation between 
K (t) and (A), we need to start with the heat kernel in a multiply-connected region. For 
two-dimensional heat kernels, there are two known results: in the simply-connected case, 
ref. provides a series expansion of the heat kernel; in the multiply-connected case, Kac 
makes a hypothesis on the heat kernel In the following, based on the Gauss-Bonnet 
theorem, we first generalize the simply-connected result of heat kernel K (t) given in Q to 
the multiply-connected case. The result can be viewed as a proof of Kac's hypothesis. 

In the simply-connected case, the heat kernel for a two-dimensional plane bounded by 
a smooth curve Tg and with the Dirichlet boundary condition on reads Q 

where S is the area of the region, L is the length of F^, k (s) is the curvature of the curve 
Tg at the point s. 

To generalize the result of the simply-connected region, eq. ([7^), to the multiply- 
connected case, we first convert the multiply-connected region bounded by F^ (figure la) 
to a simply-connected one bounded by a piecewise smooth simple closed curve F^ (figure 
lb). In the simply-connected region illustrated in figure lb, the first two terms which are 
proportional to the area and the perimeter of the region, respectively, can be calculated 
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by the method in the simply-connected case directly. For the third term, by the Gauss- 
Bonnet theorem, we can calculate the integral of the curvature along F^. The Gauss-Bonnet 
theorem reads |39] 

^{ir-ai) + j k{s)ds + jj Kda = 27rx, (7.7) 

where is the interior angle of F^ at each vertex, = in the present case is the Gauss 
curvature, and x is the Euler-Poincare characteristic number. From figure 1, since the 
integrals along 1^^ and (i = 1, • • • ,r) cancel each other in pairs, we achieve 

k{s)ds = / k{s)ds, (7.8) 

where r is the number of holes in the region. Moreover, it is also easy to see that 

r 4 

1=1 a=l 

The region bounded by F^ in figure lb is simply-connected, so x = 1- Thus, from eq. ([7.7]), 
we have 

/ k{s)ds = 2'K{l-r). (7.10) 

By the approach given by |2^, we can calculate the first two terms. Then, with eqs. ([7.6|), 
(|7!1 ), and ([7lo| ), we have 

K(t) = _L__^ + l^. (7.11) 

This is just the result that Kac hypothesized in 

By the relation between (A) and K (t), we can directly calculate the counting func- 
tion in the two-dimensional multiply-connected region: 

iV(A) = ^A--^+i^. (7.12) 
47r vr b 

8. Calculating spectrum from N (A) 

From the counting function A^(A), we can directly obtain the asymptotic expression for 
the spectrum of the system [^]. In view of the fact that one can obtain the heat kernel 
K (t) and, accordingly, obtain the counting function N (A) without knowing the spectrum 
in advance |Q, this result can serve as an approach for calculating the spectrum. 

The counting function is the number of eigenvalues smaller than A, so for the n-th 
eigenvalue A„, we have 

N{Xn) = n. (8.1) 



Therefore, from eqs. ( ^.1| ) and (4.8), we have 

fJT~+ D/2- k)'"" 



oo „ 

(47r)-^/2 y ^ ^Af/2-'= = n. (8.2) 
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The eigenvalue A„ can be solved from eq. (| 

For approximately solving eq. (^^, we assume that A„ can be expanded as 

oo 

Xn=J2a^n('--y^. (8.3) 

m=0 

Then we have 



n° + • • • . (8.4) 



4r2 (f + i) i^o^ Bo 



As examples, we list some spectra obtained by this approach for various dimensional 
balls with Dirichlet, Neumann, and Robin boundary conditions, respectively. The heat 
kernel coefficients are given in Q. 

For a three-dimensional ball, the spectrum for the Dirichlet and the Neumann or Robin 
boundary conditions is 



2^ ^ B? 8^ ' V64 J 



(8.5) 



where R is the radius. In this equation and following, the upper sign stands for the Dirichlet 
boundary condition and the lower sign for the Neumann or Robin boundary condition. 
For a four-dimensional ball, the spectrum is 



(8.6) 



(8.7) 



X - A 1/2 , 1/4 

and for a five-dimensional ball, the spectrum is 

1/ /- \2/5 1 ISvr , 1 /1125 , 20\ 1 

2 V ; i?2 32 ^ ^ R^ V 1024 3 y i?2 

9. The state density 

In this section, we provide two approaches to achieve the state density based on the above 
results. 

9.1 Calculating the state density from A^(A) and K (t) 

The first approach is straightforward. The meaning of the counting function N (A) is the 
number of the states whose eigenvalues are smaller than A, so the state density reads 
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Thus, from eq. (|4.8| ), we have 

oo „ oo 

p(A) = (4vr)-^/^ Y: YiWhk)^'""~'~" ^ ^ (4vr)-^/^5i^^/,+,5('+i)(A). 

fc=0,i,l,- '=0,1,2,- 

= (47r)-^/2 y-' ^fc A^/2-^-i + (47r)-^/2 ^fc y^D/2-k~i 

^ ' ^ T(D/2-k) ^ r(L)/2-A:) 



+ ^ (4^)-^/2i?fc5('=-f)(A) (9.2) 



From this result, one can obtain the state density once he knows the heat kernel coefficients. 
9.2 Calculating the state density from the spectrum 

An alternative way for obtaining the state density is to start with the expression of the 
counting function given in section |2|. 

Derivating both sides of eq. (|2.l| ) by A gives 



p{X) = — -f^ = lim > — ^. 9.3) 

'^^ ^ dX /3-»oo^ [e/^(A„-A) _^ lj2 ^ ^ 

Or, approximately, we can convert the summation over n to an integral: 

p(A)=lim dn — ^. (9.4) 

[e/3(A„-A) + i]2 

Then, we can directly calculate the state density from the spectrum. 



10. Conclusions 

In this paper, we reveal a relation between the counting function A'^ (A) and the heat kernel 
K (t) and provide an operable expression for N{X). By the relation, one can calculate 
the counting function from a known heat kernel, and vice versa. By the expression of the 
counting function presented in this paper, one can achieve the counting function by a direct 
calculation. 

The relation between N (A) and K (t) is an integral transformation and its inverse 
transformation. When calculating N (A) from K (t) by this relation, however, one may 
encounter the problem of divergence. This is because in most cases the expression of K (t) 
is in the form of a power series (in fact, often only first several heat kernel coefficients 
can be obtained), but the series is not uniformly convergent. As a result, when applying 
the integral transformation to the series, one needs to integrate term by term. However, 
to integrate term by term is not feasible: the integration of some terms will diverge. For 
removing the divergences, we develop a renormalization procedure. 

The results of this paper provide two approaches for calculating the counting functions: 
one is to calculate N (A) directly from the expression of N (A) given in section ^ and the 
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other is based on the relation between N (A) and K (t) given in section |3|. As a comparison 
between the two approaches, we calculate the counting function for D-dimensional boxes 
by these two approaches, respectively. 

As applications of the relation between (A) and K {t), we also calculate the count- 
ing functions for three-dimensional balls and two-dimensional multiply-connected irregular 
regions. For calculating the counting function of a two-dimensional multiply-connected 
region, we generalize the result of the heat kernel of a two-dimensional simply-connected 
region to a multiply-connected one, and the result proves Kac's conjecture on the two- 
dimensional multiply-connected heat kernel. 

Moreover, we also present approaches for calculating eigenvalue spectra and state den- 
sities from the counting function. 
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